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Abstract — Frictionless atom cooling in harmonic traps is for- 
mulated as a time-optimal control problem and a synthesis of 
optimal controlled trajectories is obtained. This work has already 
been used to determine the minimum time for transition between 
two thermal states and to show the emergence of the third law 
of classical thermodynamics from quantum thermodynamics. It 
can also And application in the fast adiabatic-like expansion of 
Bose-Einstein condensates, with possible applications in atom 
interferometry. This paper is based on our recently published 
article in SIAM J. Control Optim. Q). 

Index Terms — Quantum control, time-optimal control, atom 
cooling, quantum thermodynamics 



I. INTRODUCTION 

During the last decades, a wealth of analytical and numeri- 
cal tools from control theory and optimization have been suc- 
cessfully employed to analyze and control the performance of 
quantum mechanical systems, advancing quantum technology 
in areas as diverse as physical chemistry, metrology, and quan- 
tum information processing [2|. Although measurement-based 
feedback control (3] and the promising coherent feedback 
control [4 1 have gained considerable attention, open-loop con- 
trol has been proven quite effective. Analytical solutions for 
optimal control problems defined on low-dimensional quantum 
systems have been derived, leading to novel pulse sequences 
with unexpected gains compared with those traditionally used 
Bl- lfTOl . while numerical optimization methods, based on 
gradient algorithms or direct approaches, have been used to 
address more complex tasks and to minimize the effect of the 
ubiquitous experimental imperfections [11 1 — [ 1 3 1 . 

At the heart of modem quantum technology lies the efficient 
cooling of trapped atoms, since it has created the ultimate 
physical systems thus far for precision spectroscopy, frequency 
standards, and even tests of fundamental physics iPPfl . as 
well as candidate systems for quantum information processing 
|15|. In the present article we study a time-optimal control 
problem related to the frictionless cooling of atoms trapped 
in a time-dependent harmonic potential. Frictionless atom 
cooling in a harmonic trapping potential is defined as the 
problem of changing the harmonic frequency of the trap to 
some lower final value, while keeping the populations of 
the initial and final levels invariant, thus without generating 
friction and heating. Conventionally, an adiabatic process is 
used where the frequency is changed slowly and the system 
follows the instantaneous eigenvalues and eigenstates of the 
time-dependent Hamiltonian. The drawback of this method is 
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the long necessary times which may render it impractical. A 
way to bypass this problem is to use the theory of the time- 
dependent quantum harmonic oscillator |[T6l to prepare the 
same final states and energies as with the adiabatic process at 
a given final time, without necessarily following the instan- 
taneous eigenstates at each moment. Achieving this goal in 
minimum time has many important potential applications. For 
example, it can be used to reach extremely low temperatures 
inaccessible by standard cooling techniques [ 17] and to reduce 
the velocity dispersion and collisional shifts for spectroscopy 
and atomic clocks [ 18 1. It is also closely related to the problem 
of moving in minimum time a system between two thermal 
states |fl9l . 

It was initially proved that minimum transfer time for the 
aforementioned problem can be achieved with "bang-bang" 
real frequency controls [19]. Later, it was shown that when the 
restriction for real frequencies is relaxed, allowing the trap to 
become an expulsive parabolic potential at some time intervals, 
shorter transfer times can be obtained, leading to a "shortcut 
to adiabaticity" |20l . In our recent work [21], we formulated 
frictionless atom cooling as a minimum-time optimal control 
problem, permitting the frequency to take real and imaginary 
values in specified ranges. We showed that the optimal solution 
has again a "bang-bang" form and used this fact to obtain 
estimates of the minimum transfer times for various numbers 
of switchings. In the present article we complete our previous 
work by fully solving the corresponding time-optimal control 
problem and obtaining the optimal synthesis. The results pre- 
sented here have already been used to determine the minimum 
time for transition between two thermal states [22] and to show 
the emergence of the third law of classical thermodynamics 
from quantum thermodynamics [23], as highlighted in the 
conclusion. 

II. FORMULATION OF THE PROBLEM IN TERMS 
OF OPTIMAL CONTROL 

The evolution of the wavefunction ip(t,x) of a particle 
in a one-dimensional parabolic trapping potential with time- 
varying frequency w(i) is given by the Schrodinger equation 



ih 



~dt 



tf_d?_ ma; 2 ft) 2 
'2^8^ + 2 X 



</>, 



(1) 



where m is the particle mass and h is Planck's constant; 
i£R and ip is a square-integrable function on the real line. 
When w(t) is constant, the above equation can be solved 
by separation of variables and the solution is ip(t, x) = 
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^n( x ) are tr, e eigenfunctions of the quantum harmonic oscil- 
lator [24]. The coefficients c„ can be found from the initial 
condition c„ = f™ ^p(0,x)^f l (x)dx. 
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Fig. 1 . Time evolution of the harmonic trap frequency. 

Consider now the case shown in Fig. [T] where w(t) = ljq 
for t < and ui(t) — lot < ^o f° r t > T, This corresponds 
to a temperature reduction by a factor lut/^o, if the initial 
and final states are canonical [20|. For frictionless cooling, 
the path oj(t) between these two values should be chosen so 
that the populations of all the oscillator levels n = 0, 1, 2, . . . 
for t > T are equal to the ones at t = 0. In other words, if 



i/>(0,x) 



»=o 



Cn(0)^(x), 



and 



ip(t,x) 



^c n (i)*r(x),i>T, 

n=0 



(2) 



then frictionless cooling is achieved when 

|c„(i)| 2 Hc n (0)| 2 , i > T,n = 0, 1, 2, . . . 

Among all the paths Lv(t) that result in (fJJ, we would like 
to find one that achieves frictionless cooling in minimum 
time T. In the following we provide a sufficient condition 
on ui(t) for frictionless cooling and we use it to formulate the 
corresponding time-optimal control problem. 

Proposition 1: [T) If u/(t), with w(0) = ljq and ui(t) = 
ui(T) = lot for t > T is such that the Ermakov equation 



problem 1: Find — U\ < u(t) < u 2 with u(0) — l,u(T) = 
I/7 4 such that starting from (xi(0), #2(0)) = (1, 0), the above 
system reaches the final point {x\ (T),x 2 (T)) = (7,0),7> 1, 
in minimum time T. 

The boundary conditions on the state variables (xi, £2) are 
equivalent to those for b, 6, while the boundary conditions 
on the control variable u are equivalent to those for ui, so 
the requirements of Proposition [T] are satisfied. Parameters 
ui,U2 > define the allowable values of u(t) and it is 
u 2 > u(0) = 1. Note that the possibility w 2 (t) < (expulsive 
parabolic potential) for some time intervals is permitted [20|. 
It is natural to consider that also u\ > 1, i.e. we can at least 
achieve the negative potential V(x) — —m(jj 2 x 2 /2. Finally 
observe that the above system describes the one-dimensional 
Newtonian motion of a unit-mass particle, with position coor- 
dinate Xi and velocity x 2 . The acceleration (force) acting on 
the particle is — ux\ + l/a;f . This point of view can provide 
useful intuition about the time-optimal solution, as we will see 
later. 

In the next section we solve the following optimal control 
problem 

problem 2: Find —u\ < u(t) < u 2 , with U\,U2 > 1, such 
that starting from (xi(0),X2(0)) = (1,0), the system above 
reaches the final point (xi(T),X2(T)) — (7,0), 7 > 1, in 
minimum time T. 

In both problems the class of admissible controls formally 
are Lebesgue measurable functions that take values in the 
control set [— u\, 1*2] almost everywhere. However, as we shall 
see, optimal controls are piecewise continuous, in fact bang- 
bang. The optimal control found for problem|2]is also optimal 
for problem Q] with the addition of instantaneous jumps at 
the initial and final points, so that the boundary conditions 
u(0) = 1 and u(T) = I/7 4 are satisfied. Note that in 
connection with Fig. Q] a natural way to think about these 
conditions is that u(t) = 1 for t < and u(t) = I/7 4 for 
t > T; in the interval (0, T) we pick the control that achieves 
the desired transfer in minimum time. 



b(t)+L0 2 (t)b(t) = 



b 3 (t) 
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has a solution b(t) with 6(0) = 1, 6(0) = and b(t) = b(T) = 
(wo/wt) 1 / 2 ,^ > T, then condition (f2]i for frictionless cooling 
is satisfied. 
If we set 



X\ = 6, x 2 



b . . Lj 2 (t) 

— , U{t) = — ¥-: 



(4) 



and rescale time according to <new = ^^ol^ we obtain 
the following system of first order differential equations, 
equivalent to the Ermakov equation 



x\ = x 2l 



x 2 



—ux\ H j. 

x'l 



(5) 
(6) 



By incorporating the boundary conditions and possible restric- 
tions on io(t) due, for example, to experimental limitations, 
and setting 7 = (wo/wr) 1 ' 2 > 1, we obtain the following 
time-optimal problem for frictionless cooling 



III. OPTIMAL SOLUTION 

The system described by (O, © can be expressed in 
compact form as 

* = f(x) +ug(x), (7) 

where the vector fields are given by 
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(8) 



and x G V = {(xi,x 2 ) G R 2 : Xi > 0} and u G 
U = [—ui,u 2 ]. Admissible controls are Lebesgue measur- 
able functions that take values in the control set U. Given 
an admissible control u defined over an interval [0,T], the 
solution x of the system ((7j corresponding to the control u is 
called the corresponding trajectory and we call the pair (x, u) 
a controlled trajectory. Note that the domain T> is invariant in 
the sense that trajectories cannot leave T>. Starting with any 
positive initial condition xi(0) > 0, and using any admissible 
control u, as x\ —¥ + , the "repulsive force" l/xf leads to 



an increase in x\ that will keep x\ positive (as long as the 
solutions exist). 

For a constant Ao and a row vector A = (Ai, A2) € (R 2 ) 
define the control Hamiltonian as 

H = H(X ,X,x,u) = A + (X,f(x) +ug(x)). 



Pontryagin's Maximum Principle [25 1 provides the following 
necessary conditions for optimality: 

Theorem 1 (Maximum principle): 11251 Let (x*(t),u*(t)) 
be a time-optimal controlled trajectory that transfers the initial 
condition x(0) = xq into the terminal state x(T) = xt- Then 
it is a necessary condition for optimality that there exists 
a constant Ao < and nonzero, absolutely continuous row 
vector function X(t) such that: 

1) A satisfies the so-called adjoint equation 

A(t) = - — (Ao,A(*),x.(*),u*(t)) 

ox 

2) For < t < T the function u h->- H(Xo, \(t),x*(t),u) 
attains its maximum over the control set U at u = u*(t). 

3) jy(Ao,A(t),a?*(t),u.(t))sO. 

We call a controlled trajectory (x, u) for which there exist 
multipliers Ao and X(t) such that these conditions are satisfied 
an extremal. Extremals for which Ao = are called abnormal. 
If Ao < 0, then without loss of generality we may rescale the 
A's and set Ao = — 1. Such an extremal is called normal. 

For the system (0, © we have 



if (A , A, x, u) = A + X\x 2 + A 2 



and thus 



A = -A 







1 



-(u + 3/xf) 



X\U 



= -XA 



(9) 



(10) 



Observe that H is a linear function of the bounded control 
variable u. The coefficient at u in H is — X 2 x\ and, since 
Xi > 0, its sign is determined by $ = — A2, the so-called 
switching function. According to the maximum principle, point 
2 above, the optimal control is given by u = — u% if $ < and 
by u — u 2 if $ > 0. The maximum principle provides a priori 
no information about the control at times t when the switching 
function $ vanishes. However, if $(i) =0 and $(t) ^ 0, then 
at time t the control switches between its boundary values and 
we call this a bang-bang switch. If $ were to vanish identically 
over some open time interval I the corresponding control is 
called singular. 

Proposition 2: For Problem [2] optimal controls are bang- 
bang. 

Proof: Whenever the switching function $(£) = — X2(t) 
vanishes at some time t, then it follows from the non-triviality 
of the multiplier X(t) that its derivative $(£) = — X2(t) = 
Xi(t) is non-zero. Hence the switching function changes sign 
and there is a bang-bang switch at time t. ■ 

Definition 1: We denote the vector fields corresponding to 
the constant bang controls — m and u 2 by X = f — U\g 
and Y = f + u 2 g, respectively, and call the trajectories 
corresponding to the constant controls u = —U\ and u = u 2 
X- and F-trajectories. A concatenation of an X-trajectory 



followed by a F-trajectory is denoted by XY while the 
concatenation in the inverse order is denoted by YX. 

In this paper we establish the precise concatenation se- 
quences for optimal controls and in particular calculate the 
times between switchings explicitly. 

Proposition 3: All the extremals are normal. 

Proof: If (x, u) is an abnormal extremal trajectory that 
has a switching at time t, then, since X 2 (t) — 0, it follows 
from H — that we must have x 2 (t) = 0. The starting point 
is (1,0) and suppose that u = —u\ initially. From © it is 
x 2 > so x 2 > and a switching at a point with x 2 (t) > 0, 
not allowed for an abnormal extremal, is necessary in order to 
reach the target point (7, 0). If u = u 2 initially, then ±2(0) = 
1 — u 2 < and x 2 < for some time interval. During this 
time it is xi < and consequently x\ < 1 < 7. A switching 
is necessary, which takes place on the xi-axis for an abnormal 
extremal. The control changes to u = —u\ and the situation 
is as before, where one more switching is necessary at a point 
with x 2 (t) > 0, forbidden for abnormal extremals. Thus, there 
are no abnormal extremals in the optimal solutions. ■ 

For normal extremals we can set A = — 1. Then, H = im- 
plies that for any switching time t we must have Xi(t)x 2 (t) = 
1. For an XY junction we have 4>(t) = X\(t) > and thus 
necessarily x 2 (t) > and analogously optimal YX junctions 
need to lie in {x 2 < 0}. We now develop the precise structure 
of the switchings in a series of Lemmas. 

Lemma 1 (First integrals): A first integral of the motion 
along the X-trajectory passing through (a, 0) is 



2 1 2 1 

uix-, H — 7j = -ma H — o, 
x\ 



■n 



a- 



(11) 



while a first integral of the motion along the Y- trajectory 
passing through (/3, 0) is 

1 .0 1 



u 2 x\ 



u 2 p' 



(12) 



*i P 

Proof: Use the system equations (0 and ©. ■ 

Lemma 2 (Inter-switching time): Let p = (xi,x 2 ) be a 
switching point and r denote the time to reach the next 
switching point q. If pq is a y-trajectory, then 



sin(2^/u2"T) 



2Ju^x\X 2 



x 2 + u 2 x\ 

while, if pq is an X -trajectory, then 

2Ju{xix 2 



, cos(2 v /u 2 "r) 



u 2 xl 



u 2 x 



(fe) 



r . cosh(2 v / uir) 

U\X\ 



U\X\ 



U\X 



sinh(2y / U] _ T) = — - 

"" "" "1^1 

(14) 

Note that the inter-switching times depend only on the ratio 

x 2 /x 1 . 

Proof: These formulas are obtained as an application of 

the concept of a "conjugate point" for bang-bang controls 

|26|. Without loss of generality assume that the trajectory 

passes through p at time and is at q at time r. Since 

p and q are switching points, the corresponding multipliers 

vanish against the control vector field g at those points, i.e., 

(X(0),g(p)) = (X(T),g(q)) = 0. We need to compute what 

the relation (X(r),g(q)) = implies at time 0. In order to do 

so, we move the vector g(q) along the y-trajectory backward 



from q to p. This is done by means of the solution w(t) of 
the variational equation along the y-trajectory with terminal 
condition w(t) = g(q) at time r. Recall that the variational 
equation along Y is the linear system w = Aw where A 
is given in ( TTOb . Symbolically, if we denote by e tY (p) the 
value of the ^-trajectory at time t that starts at the point p 
at time and by (e~* Y )» the backward evolution under the 
linear differential equation w = Aw, then we can represent 
this solution in the form 

«;(0) = (e- rY )Mr) = (e- TY ),g(q) 

= (e- TY U(e TY (p)) = (e~ TY )* o g o e^ Y (p). 

Since the "adjoint equation" of the Maximum Principle is 
precisely the adjoint equation to the variational equation, it 
follows that the function t h-> (X(i),w(i)) is constant along 
the ^-trajectory. Hence (A(r),g(q)) = implies that 

(A(0), W (0)) = (\(0),(e- TY U(e TY (p)))=Q 

as well. But the non-zero multiplier A(0) can only be or- 
thogonal to both g(p) and w(0) if these vectors are parallel, 
g(p)\\w(0) — {e~ TY )*g{e TY (p)). It is this relation that defines 
the switching time. 

It remains to compute w(0). For this we make use of the 
well-known relation ll27l 



(e 



- TY ).°go? Y 



idY 



(.9) 



(15) 



where the operator adY is defined as adY(g) = [Y,g], with 
[, ] denoting the Lie bracket of the vector fields Y and g. For 
our system, the Lie algebra C generated by the fields / and g 
actually is finite dimensional: we have 



[/,<?](*) 



-xi 



and the relations 

[/,[/,<?]] =2/, [g,[f,g}} = -2g 

can be directly verified. Using these relations and the analyt- 
icity of the system, e tadY (g) can be calculated in closed form 
from the expansion 



e ta dY{g) = J2 d n Y{g) 



(16) 



n=0 



where, inductively, ad n Y(g) = [Y,ad n 1 Y(g)]. It is not hard 
to show that for n = 0, 1, 2, . . ., we have that 



and 



so that 



ad 2n+1 Y(g) = (-4u 2 T[f,g] 



ad^Y(g) = 2(-4u 2 ) n (f-u 2 g), 
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n=0 v '' 
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Fig. 2. Consecutive switching points lie on two opposite-slope lines through 
the origin. Blue dashed curves correspond to X-segments, red dotted curves 
to Y-segments. 



By summing the series appropriately we obtain 

1 



tadY 



(g) =g + 



2JT^ 



x(2^t)[f,g] 



1 

2u : 



[1 - cos(2y/itot)](f - u 2 g). 



Hence the field w(0) — (e TY )*g(e rY (p)) is parallel to 
g(p) = (0, — xi) T if and only if 



u 2 xi sin(2v / ujr) + x 2 [1 — cos(2 y /w 2 T)] = 0. 



Hence 



sin(2 v A uJr) 



x 2 



U 2 Xi 



■[1 - cos(2\/w^t)] 



(17) 



from which ( fT~3T > follows. Note that the solution cos(2^/u 2 t) = 
1 is rejected because it corresponds to r = or r = tt/^/u^, 
the latter being the period of the closed trajectory. 

Working analogously for an X-trajectory we obtain ([T4V ■ 
Lemma 3 (Main technical point): The ratio of the coordi- 
nates of consecutive switching points has constant magnitude 
but alternating sign, while these points are not symmetric with 
respect to the xi-axis. 

Proof: Consider the trajectory shown in Fig. [2] with 
switching points (k,/i), (C, £) and (A, v). Starting from (k,//) 
and integrating the equations of motion <(5j and © for the 
inter-switching time given in < fT~3T >, we can find the coordinates 
of the next switching point and then show that £/£ = — /x/« 
while ((,£,) y^ (k, — /x) JT]. Subsequently, integrating the 
equations for the inter-switching time given in (fl4t . we can 
also show that vj\ = -£/( and (A, v) ^ (C, -£) Q. 

Alternatively, we present a more elegant proof based on 
the symmetries of the system. Observe that the transformation 
(t,Xi,x 2 ) —> (—t,xi,—x 2 ) leaves the system © and © 
invariant for constant u. So, starting from (£, ^) and running 
the system backwards to the next (previous in the forward 
direction) switching point (k, /j), the switching time is given 
by a relation similar to ([PJl l. with x 2 jx\ — —£/C B ut this 
switching time is the same as in the forward direction, where 
x 2 /x\ — h/k in (TT~3T >- Then, using (fT3] l, it is not hard to see 
that £/£ = — h/k. Note that C = k would imply ^ = — /j,, 
i.e. returning to the same point on the iri-axis with opposite 
velocity before switching again, which is obviously not time 
optimal. Thus, it is £ 7^ n in general so (£,£) 7^ (k, — ju). 
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Fig. 3. Blue dashed curves correspond to X-segments, red dotted curves to 
y-segments. 



In the following proposition we use Lemma [3] to determine 
the form of the optimal trajectory. 

Proposition 4 (Form of the optimal trajectory): The opti- 
mal trajectory can have the one-switching form XY or the 
spiral form YX . . . YXY with an even number of switchings. 
Proof: We first show that when the optimal trajectory has 
more than one switching, it cannot start with an X-segment. 
For just two switchings, consider the trajectory XYX depicted 
in Fig. |2 where a = 1 (starting point), (7,0), 7 > 1 is the 
target point and (k,/x), (Cj£) are tne switching points. Since 
both of the switching points belong to the F-segment passing 
through (f3, 0), their coordinates satisfy (fT2t . If we denote by 
s the common ratio p 2 /k 2 = £ 2 /C 2 = s, then both k, £ satisfy 
the equation 



(a) 



{s + u 2 )x{ - [u 2 p 2 + — )x\ 



1 = 0, 



so 



K 2 C 2 = 



<1, 



U2 

2a2 



since u 2 > l,s > 0. But also k £ > 1, since k > 1 and 
C 2 > 7 2 > 1. Thus this trajectory cannot be optimal. 

For more switchings, consider the case shown in Fig |2] 
where now a = 1, and use s to denote the common ratio of 
the squares of the coordinates at the switching points. If r is 
the switching time between ((, £) and (A, v), then from (fl4] i 

we obtain 

s cosh(2rY / u 1 ) + 1 



> 1. 



< 1 <=> (tii - 1 



)K 2 > 



1, 



cosh(2r v /u7) — 1 
But s/ui = fi 2 /(uiK 2 ), and from (fTTT i we find (a = 1) 

s (um 2 + 1)(k 2 - 1) 

since ui > 1. Thus if the optimal trajectory has more than 
one switching, it needs to start with a F-segment. 

We next show that the optimal trajectory reaches the target 
point (7,0), 7 > 1 with a Y"-segment. This is obviously the 
case for one switching, and also for two switchings since only 
the YXY trajectory is permitted (the XYX was excluded 
above). For more than two switchings consider the situation 
shown in Fig. [3] It is /1 2 / 'k 2 = £ 2 /C, 2 = s and s > u\ since 
at least one Y"XY"-segment is included in the trajectory. Point 
(£, £) belongs to the final X-segment ending to (7, 0), so 

1 „ 1 



(s-u 1 )( 2 + 



c 2 



-«i7 



\ (ft 



Fig. 4. 
turns. 



(b) 
(a) Trajectory with one switching (zero turns) (b) Trajectory with n 



The left hand side is positive, since s > u\, while the right 
had side is negative, since 7 > 1,U\ > 1. Thus the optimal 
trajectory reaches the target point with a Y"-segment. ■ 

Corollary 1: For \u\ < 1 the optimal solution has only one 
switching. 

Proof: For u = u 2 = 1 the starting point (1,0) is an 
equilibrium point of system (0, ©. So the optimal trajectory 
cannot start with a Y"-segment. The only trajectory thus 
permitted is XY ■ 

From Proposition |4] we see that the optimal trajectory can 
have aside from the expected one-switching form, shown in 
Fig. |4(a)| the spiral form shown in Fig. |4(b)| An intuitive 
understanding of this latter form can be obtained by viewing 
system equations (0, © as describing the motion of a unit 
mass particle with position X\ and velocity x 2 . In light of 
this interpretation we see that along a spiral trajectory the 
particle, instead of moving directly to the target, goes close to 
x\ = where there is a strong repulsive potential (1/xf) to 
acquire speed and reach the target point faster. In the following 
theorem we calculate the transfer time for the candidate 
optimal trajectories. 

Theorem 2: Starting from (1,0), the necessary time to 
reach the target point (7,0), 7 > 1 with one switching is 



T = 



1 



sinh 



/ mi(7 2 - 1)(m 2 7 2 - 1) 
7 2 (ui +u 2 )(u 1 + 1) ^ 

f u 2 (7 2 -l)(m 7 2 + l) \ 

(ui +u 2 ){u 2l i - 1) I 



(18) 



The necessary time to reach the target with n turns (2n 
switchings) is 



r 



T n = T! + nT x + {n- l)T Y + T F , 



(19) 



where 



T, 



1 



T F 



1 
2Jui 



: COS 



sci +u 2 \Jc\ -4(. 


H-tfcM 


(S + M 2 )v /C l- 


1 — 
4u 2 / 


-SC„ + l+M2yC ? 2 l+ l 


-4(s + u 2 ) 



(20) 



(s + -u 2 )y / 4+i - 4u 2 



Tx = 



1 



Ty = 



2yui 

1 

2^/M 2 



: COsh 



2tt 



cos 



S + Ml 

s — U\ 
, I s-u 2 



u 2 



c\ = u 2 + 1, 



1 



u 2 7 + -o, 
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and s is the solution of the transcendental equation 

s — u\ 



ci + vq " 4 ( s + "2) 
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(24) 
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in the interval ui < s < (w 2 — l) 2 /4. The constants 
c\ and c„+i characterize the first and the last F-segments, 
respectively, of the trajectory. The number of turns satisfies 
the following inequality 

To 



n < 



(27) 



Tx(s+)_ 

where s + = (u 2 — l) 2 /4 and [] denotes the integer part. 

Proof: In Fig. |4(a)| we show a trajectory with one 
switching point B(k, /i). The coordinates of this point satisfy 
equations ( fTTT i and (fT2l with a = 1 and f3 = 7, 
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1 9 1 



fll~UlK 2 -\ - = "l-Ul, 



fl + U 2 K 

from which we find 



w 2 7 
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u 27 4 + l+7 2 (ui-l) 



7 2 (ui + u 2 ) 

By integrating the equations of motion (0 and ©, we find 
the necessary time along each segment of the trajectory, AB 
and BC. The total transfer time is the sum of these times 
and is given by (1 1 8t> . Next consider the case with n turns and 
2n switching points (kj,Hj), Fig. |4(b)| with constant ratio 
/i 2 //« 2 = s. The first switching point satisfies the equations 

-^ = ci, (28) 

(29) 



M? 



u 2 n\ 
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,u 2 - UlK 2 H ^ = c, 

«1 



where ci is given by (124b and c 
second switching point satisfies 



-u\a\ 



A 



A 



u 2 k 2 + — = c 2 , 

K 2 



UiK 2 



c, 



1/af, while the 

(30) 
(31) 



where c 2 = u 2 /3f + l/Pf. The constants c\ and c 2 characterize 
the first and second F-segments of the trajectory, while the 
constant c characterizes the X-segment joining them. Sub- 
tracting (|29l from (|3H and using Lemma [3] which assures 
that Ki y^ k 2 (consecutive switching points are not symmetric 
with respect to xi-axis) we find that 

ui ^-^j = 0. (32) 



/v-i f\"-\ 



But from d28l . (T30b and the constant ratio relation we find 

2 



c i + V c l - 4 ( s + u 2) 
2 _ c 2 + y/^ - 4(s + u 2 ) 

2(s + u 2 ) 

where, while solving the quadratic equations we used the — 
sign for the first and the + sign for the second switching 
point. The choice of sign for the first switching point will 
be justified below, while the choice of sign for consecutive 
switching points should be alternating to avoid picking the 
symmetric image of the previous point. Using these relations, 
takes the form 



c\ 



+ Vq 1 - 4(s + u 2 ) 



s — U\ 



C 2 + y/(% - 4(s + u 2 ) s + u 2 

By repeating the above procedure for all the consecutive pairs 
of switching points, we find 



+ y/<% - 4(a + us) 



111 



Ci+l 



'i+l 



4(s + u 2 ) 



u 2 



i = 1,2, . . . ,n. 



Multiplying the above equations we obtain d26l l, one transcen- 
dental equation for the ratio s. If we choose the + sign in 
the quadratic equation for k\, we obtain an equation similar 
to d26b but with inverted left hand side. It is c n+ i > c\ <=> 

(7 2 - l)(u 2 7 2 - 1) > and c 1: c n+1 > 0, so 



n+l + J C 2 n+1 - A{S + U 2 ) 



Cl + y/4-Ms + U2) 



> 1 > 



"1 



U 2 



and the corresponding transcendental equation has no solution. 
Note that the left hand side of d26i > is a decreasing function 
of s while the right hand side is an increasing one, so if a 
solution exists, it is unique. The ratio is bounded below by 
the requirement s/u\ > 1 and above by c 2 — 4(s + u 2 ) > 
<=> s < s+ = (u 2 — l) 2 /4. This is also the maximum value 
of s on the first F-segment ( f28l . Once we have calculated 
this ratio, we can find the time interval between consecutive 
switchings using ( f22l for an X-segment and ( f23l for a Y- 
segment, relations obtained from Lemma |2] on the inter- 
switching time. Observe that the times along all intermediate 
X- (respectively Y-) trajectories are equal. The initial time 
interval Tj (from the starting point up to the first switching) 
and the final time interval Tp (from the last switching up 
to the target point) can be easily calculated and are given 
in (l20b and (|2TT >. respectively. The total duration T n of the 
trajectory with n turns joining the points (1,0) and (7,0) is 
given by (QjJJ. Observe that T n (s) > nT x (s) > nT x (s+), 
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Fig. 5. (a) Transfer times corresponding to zero, one and two turns for Fi §- 6 - ( a ) Transfer times corresponding to zero, one, two and three turns 
U1 = 1,« 2 = 8,7 G [1,10]. (b) Switching curves (black solid curves) and io \ u i = }'^ = 50 ' 1 G [ l . ^.(b) Switching curves (black solid curves) 
characteristic optimal trajectories starting from (1,0). 



and characteristic optimal trajectories starting from (1,0). 



where the last inequality follows from the fact that Tx is a 
decreasing function of s, see (1221 1, A solution with n turns 
can be candidate for optimality only if the number of turns is 
bounded as in d27| i. Otherwise we have T n (s) > Tq and the 
one-switching strategy is faster. ■ 

Using Theorem [2] we can find the times T n for a specific 
target (7, 0) and compare them to obtain the minimum time. 

IV. EXAMPLES 



In Fig. |5(a)| we plot the times Tq, T\ and T-i from Theorem|2] 
corresponding to zero, one and two turns, for u\ = l,ti2 — 8 
and 7 e [1, 10]. For 7 < 71 the strategy with zero turns (one 
switching) is optimal, while for 7 > 71 it is the strategy with 
one turn (up to the range of 7 plotted). The point (71, 0) can be 
reached with both strategies in equal time, that is, it belongs 
to the cut-locus of these two control sequences from (1,0). 
Note that the strategies with one and two turns are feasible 
after some 7 > 1, where the transcendental equation (|26| | has 
a solution. In Fig. |5(b)| we plot the switching curves (black 
solid curves) as well as some characteristic optimal trajectories 
starting from (1,0). For 7 < 71 the optimal trajectory starts 
with an A-segment that coincides with the switching curve 
(black solid curve) passing from (1,0). It switches at some 
point and then travels along a F-segment (red dotted curve) 
to meet the ii-axis. For 7 > 71 the optimal trajectory starts 
with a F-segment (red dotted curve passing from (1,0)) and 
switches at some point in the tiny black area of this curve 
to an X-segment (blue dashed curve). Then it meets at some 
point the second switching curve on the upper quadrant and 
changes to a F-segment (red dotted curve) that hits the xi-axis 
at the target point. Note that the optimal trajectories between 
the two switchings (blue dashed curves) are very close to the 



second switching curve on the upper quadrant and they are 
not shown entirely. 

In Fig. |6(a)| we plot the times T ,T 1 ,T 2 and T 3 from 
Theorem |2] corresponding to zero, one, two and three turns, 
for u\ = 1,U2 — 50 and 7 € [1, 15]. Again, for small 7 the 
one-switching strategy is optimal and after some 7 = 71 the 
one-turn strategy becomes faster, but there is also a 7 = 72 
beyond which the two-turn strategy is optimal (up to the range 
of 7 plotted). The point (72,0) thus belongs to the cut-locus 
of the one- and two-turn control sequences from (1,0) since it 
can be reached with one or two turns in equal time. In Fig. |6(b)| 
we plot the switching curves (black solid curves) along with 
some characteristic optimal trajectories starting from (1,0). 
For 7 > 72 the optimal trajectory makes an additional turn. 
This is demonstrated by the three adjacent ^-segments (red 
dotted curves), which switch close to to the corresponding 
X-segments (blue dashed curves), on a tiny switching curve 
which is hardly seen. In turn, these trajectories switch on the 
third switching curve on the upper quadrant to y-segments 
(red dotted curves) that hit the rri-axis at the target points. 



V. CONCLUSION AND FUTURE WORK 

In this article we solved a time-optimal control problem 
related to frictionless atom cooling in harmonic traps. The 
results presented here can be immediately extended to the 
fast frictionless expansion of a two-dimensional Bose-Einstein 
condensate confined in a parabolic trapping potential [28 1, 
with possible application in atom interferometry, and even to 
the implementation of a quantum dynamical microscope, a 
controlled expansion that allows to scale up an initial many- 
body state of an ultracold gas by a desired factor while 
preserving the quantum correlations of the initial state ll29ll . 



This work has also been used to show how the third 
law of classical thermodynamics, known as unattainability 
principle, emerges from quantum thermodynamics [23]. In a 
dynamical interpretation, this law states that absolute zero is 
unattainable, since the cooling rate from a thermal bath with 
falling temperature declines as well and approaches zero with 
an appropriate power of the temperature. The heat machine 
used to demonstrate this is a quantum refrigerator, the quantum 
analog of the classical Otto cycle, where the working medium 
is made up of particles in a harmonic (possibly repelling) 
potential instead of classical particles in a piston. The initial 
frequency coincides with the maximum allowed frequency as 
well as the strongest repelling frequency, so u\ = 112 = 1. 
From Corollary Q] we see that the minimum time T for the 
adiabatic-like cooling branch is given by equation ( fTST l in 
Theorem 12 As the temperature approaches zero, wy — > and 
7 = y/ujo/ujr ^ > 00, so the cooling time approaches infinity 
logarithmically, T ~ In 7 |22|. 
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